Repulsive, nonmonotonic Casimir forces in a glide-symmetric geometry 
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We describe a three-dimensional geometry that exhibits a repulsive Casimir force using ordinary 
metallic materials, as computed via an exact numerical method (no uncontrolled approximations). 
The geometry consists of a zipper-like, glide- symmetric structure formed of interleaved metal brack- 
ets attached to parallel plates. Depending on the separation, the perpendicular force between the 
plates/brackets varies from attractive (large separations) to repulsive (intermediate distances) and 
back to attractive (close separations), with one point of stable equilibrium in the perpendicular 
direction. This geometry was motivated by a simple intuition of attractive interactions between 
surfaces, and so we also consider how a rough proximity-force approximation of pairwise attractions 
compares to the exact calculations. 



In this letter, we describe a metallic, glide-symmetric, 
"Casimir zipper" structure (depicted in Fig. [I]) in which 
both repulsive and attractive Casimir forces arise, includ- 
ing a point of stable equilibrium with respect to per- 
pendicular displacements. We compute the force using 
an "exact" computational method (i.e. with no uncon- 
trolled approximations, so that it yields arbitrary accu- 
racy given sufficient computational resources), and com- 
pare these results to the predictions of an ad hoc at- 
tractive interaction based on the proximity-force approx- 
imation (PFA). Casimir forces, a result of quantum vac- 
uum fluctuations, arise between uncharged objects, most 
typically as an attractive force between parallel metal 
plates pQ that has been confirmed experimentally [2j [3] . 
One interesting question has been whether the Casimir 
force can manifest itself in ways very different from this 
monotonically decaying attractive force, and especially 
under what circumstances the force can become repul- 
sive. It has been proven that the Casimir force is al- 
ways attractive in a mirror-symmetric geometry (with 
e > 1 on the imaginary- frequency axis) [4 , but there 
remains the possibility of repulsive forces in asymmetric 
structures. For example, repulsive forces arise in exotic 
asymmetric material systems, such as a combination of 
magnetic and electric materials [5j [6j [7] , fluid-separated 
dielectric plates [8], metamaterials with gain [9], or ex- 
cited atoms [10]. Another route to unusual Casimir phe- 
nomena is to use conventional materials in complex ge- 
ometries, which have been shown to enable asymmetri- 
cal lateral "ratchet" effects [11] and nonmonotonic de- 
pendencies on external parameters [12]. Until recently, 
however, predictions of Casimir forces in geometries very 
different from parallel plates have been hampered by the 
lack of theoretical tools capable of describing arbitrary 
geometries, but this difficulty has been addressed by re- 
cent numerical methods [13] [14] [15] [16]. In this letter, 
we use a technique based on the mean Maxwell stress 
tensor computed numerically via an imaginary-frequency 
Green's function, which can handle arbitrary geometries 
and materials [T3j . 

The geometry that we consider is depicted schemati- 




FIG. 1: (Color) Three-dimensional schematic of the Casimir 
"zipper" geometry of interlocking metal brackets (shown 
in different colors for illustration only), along with a two- 
dimensional xy cross-section. The dashed lines extruding 
from the plates to the squares indicate their out-of-plane con- 
nectivity. 



cally in Fig. [T] we have two periodic sequences of metal 
"brackets" attached to parallel metal plates, which are 
brought into close proximity in an interlocking "zip- 
per" fashion. In Fig. [I] we have colored the two 
plates/brackets red and blue to distinguish them, but 
they are made of the same metal material. This structure 
is not mirror symmetric (and in fact is glide-symmetric, 
although the glide symmetry is not crucial), so it is not 
required to have an attractive Casimir force by Ref. [4j 
Furthermore, the structure is connected and the ob- 
jects can be separated via a rigid motion parallel to the 
force (a consideration that excludes interlocking "hooks" 
and other geometries that trivially give repulsive forces) . 
This structure is best understood by considering its two- 
dimensional cross-section, shown in Fig. Wright) for the 
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middle of the brackets: in this cross-section, each bracket 
appears as an s x s square whose connection to the ad- 
jacent plate occurs out-of-plane. (Here, the brackets are 
repeated in each plate with period A = 2s + 2h and are 
separated from the plates by a distance d. The plates are 
separated by a distance 2d + s + a, so that a = is the 
point where the brackets are exactly aligned.) The mo- 
tivation for this geometry is an intuitive picture of the 
Casimir force as an attractive interaction between sur- 
faces. When the plates are far apart and the brackets are 
not interlocking, the force should be the ordinary attrac- 
tive one. As the plates move closer together, the force 
is initially dominated by the attractions between adja- 
cent bracket squares, and as these squares move past one 
another (a < in Fig. [I]), one might hope that their at- 
traction leads to a net repulsive force pushing the plates 
apart. Finally, as the plates move even closer together, 
the force should be dominated by the interactions be- 
tween the brackets and the opposite plate, causing the 
force to switch back to an attractive one. This intu- 
ition must be confirmed by an exact numerical calcu- 
lation, however, because actual Casimir forces are not 
two-body attractions and can sometimes exhibit qualita- 
tively different behaviors than a two-body model might 
predict [17]. Such a computation of the total force per 
unit area is shown in Fig. [2j and demonstrates precisely 
the expected sign changes in the force for the three sep- 
aration regimes. These results are discussed in greater 
detail below. 

Previous theoretical studies of Casimir forces in ge- 
ometries with strong curvature have considered a va- 
riety of objects and shapes. Forces between isolated 
spheres [16] and isolated cylinders [18 , or between a sin- 
gle sphere [19], or cylinder [TSJ [19] and a metal plate, 
all exhibit attractive forces that decrease monotonically 
with separation. When a pair of squares [12 or cylin- 
ders [18 interacts in the presence of two adjacent metal 
sidewalls, the force is still attractive and monotonic in 
the square/square or cylinder /cylinder separation, but 
is a nonmonotonic function of the sidewall separation. 
When two corrugated surfaces are brought together in a 
way that breaks mirror symmetry (i.e., the corrugations 
are not aligned between the two surfaces), a lateral force 
can arise [20j [21] , and an asymmetric lateral force from 
asymmetric corrugations can lead to a "ratchet" effect 
in which random forces preferentially displace the plates 
in one direction [11]. Such a lateral force has also been 
observed experimentally [22]. In the geometry of Fig. [I] 
in contrast, there is no lateral force (due to a mirror- 
symmetry plane perpendicular to the plates), and hence 
we consider only the normal force between the plates. 
Because of the strong curvature of the surfaces relative 
to their separations, simple parallel-plate approximations 
are not valid (although we consider their qualitative ac- 
curacy below), and the force must be computed numeri- 
cally. 
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FIG. 2: (Color) Top: Plot of the Neumann (blue, TE), Dirich- 
let (red, TM) and total (black, TE+TM) Casimir pressure (in 
units of he/ A 4 ) between the objects of Fig. [I] as a function 
of a/s. The inset illustrates a two-dimensional cross-section. 
Bottom: Schematic indicating the various qualitatively differ- 
ent Casimir force regimes between the two structures. 



The numerical method we employ is based on inte- 
gration of the mean stress tensor, evaluated in terms 
of the imaginary-frequency Green's function via the 
fluctuation- dissipation theorem [13]. The Green's func- 
tion can be evaluated by a variety of techniques, but 
here we use a simple finite-difference frequency-domain 
method [I3j [23] that has the advantage of being very 
general and simple to implement at the expense of 
computational efficiency. In particular, the computa- 
tion involves repeated evaluation of the electromagnetic 
Green's function, integrated over imaginary frequency 
w = —iuo and a surface around the object of interest. 
The Green's function is simply the inverse of a linear op- 
erator [V x V x -\-w 2 e(iw,r)], which here is discretized 
using a finite-difference Yee grid [23] and inverted us- 
ing the conjugate- gradient method [24 . In order to sim- 
plify the calculations, we assume the length of the brack- 
ets in the z direction L to be sufficiently long to make 
their contributions to the force negligible (we estimate 
the minimum length below). We can therefore describe 
the geometry as both z-invariant and ^/-periodic (with 
period A). This implies that it is only necessary to com- 
pute the Green's function using an xy unit cell, with the 
periodic /invariant directions handled by integrating over 
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the corresponding wavevectors [13]. Furthermore, we ap- 
proximate the bracket /plate materials by perfect metals, 
valid in the limit of small lengthscales (which are dom- 
inated by long- wavelength contributions where the skin 
depth is negligible). In this case, the contributions to 
the force can be separated into two polarizations: trans- 
verse electric (TE) with the electric field in the xy plane 
(a scalar magnetic field with Neumann boundary condi- 
tions); and transverse magnetic (TM) with the magnetic 
field in the xy plane (a scalar electric field with Dirichlet 
boundary conditions) [13], and these two contributions 
are shown separately in Fig. [2] 

The resulting force per unit area between the plates, 
for the chosen parameters d/s = 2 and h/s = 0.6, is 
plotted as a function of a/s in Fig. [5] (Top); error bars 
show estimates of the numerical accuracy due to the fi- 
nite spatial resolution. A number of unusual features are 
readily apparent in this plot. First, the sign of the force 
changes not only once, but twice. The corresponding ze- 
ros of the force lie at a/s « —0.8 and a/s « — 10 -2 . 
The first zero, a/s ~ —0.8, is a point of unstable equilib- 
rium, to the left of which the force is attractive and to 
the right of which the force is repulsive. The second zero 
at a/s « — 10 -2 corresponds to a point of stable equi- 
librium, with respect to perpendicular displacements, for 
which the force is attractive to the right and repulsive 
to the left. (This point is still unstable with respect to 
lateral displacements, parallel to the plates and perpen- 
dicular to the brackets, however: any such lateral dis- 
placement will lead to a lateral force that pulls the red 
and blue brackets together.) In between these equilibria, 
the repulsive force has a local maximum at a/s ~ —0.5. 
Finally, at a/s ~ 0.6 the magnitude of the attractive 
force reaches a local maximum (a local minimum in the 
negative force on the plot), and then decreases asymptot- 
ically to zero as a/s — > oo. Thus, as the two objects move 
apart from one another, the force between them varies in 
a strongly nonmonotonic fashion (distinct from the non- 
monotonic dependence on an external parameter shown 
in our previous work [12j [TTj [18] ) . These three different 
sign regimes are shown schematically in Fig. [2] (Bottom) , 
as predicted by the intuitive picture described above. 

Since the qualitative features of the Casimir force in 
this geometry correspond to the prediction of an intuitive 
model of pairwise surface attractions, it is reasonable to 
ask how such a model compares quantitatively with the 
numerical results. The most common such model is the 
proximity- force approximation (PFA), which treats the 
force as a summation of simple "parallel-plate" contri- 
butions [19] . (Another pairwise power-law heuristic is 
the "Casimir-Polder interaction" approximation, strictly 
valid only in the limit of dilute media [25 ].) Applied to 
a geometry with strong curvature and/or sharp corners 
such as this one, PFA is an uncontrolled approximation 
and its application is necessarily somewhat ad hoc (due 
to an arbitrary choice of which points on the surfaces 



to treat as "parallel plates"), but it remains a popular 
way to quantify the crude intuition of Casimir forces as 
pairwise attractions. 
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FIG. 3: (Color) Comparison of Casimir pressure (in units 
of he/ A 4 ) as a function of a/s between the stress-tensor (ex- 
act) numerical results (black squares) and the proximity-force 
approximation (solid green). Also shown are the individual 
square-square (dashed blue) and square-plate (dashed or- 
ange) contributions to the PFA force. Inset: Schematic il- 
lustration of the chosen PFA "lines of interaction" between 
squares (dashed black lines). 

Applying the PFA approximation to the two objects 
in Fig. [I] we treat the net force as a sum of three con- 
tributions: the force between the two parallel plates, the 
force between each square and the opposite plate, and the 
force between adjacent red and blue squares. Namely, 

Py>fa = (^pi-pi + 2F S q_pi + 2F sq _ sq ) , (1) 

where the first term is the pressure between two parallel 
plates (pl-pl), and the two remaining terms correspond 
to the square-plate (sq-pl) and square-square (sq-sq) 
interactions. The factors of A and L are introduced be- 
cause these expressions are computed per unit length in 
the z direction, and per period in the y direction. 

The first two PFA contributions are relatively simple 
to calculate because they are between parallel metal sur- 
faces, and thus (in the PFA approximation) are the or- 
dinary Casimir force weighted by the respective areas: 



pl " pl 120A (2d + a + s) 4 v 7 

hC7T 2 S 1 

sq - pl = " 240A Jd + aY (3) 
Computing the square-square force is less straightfor- 
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ward, since there is some ambiguity as to what the PFA 
approximation even means for two non-parallel surfaces 
(separate from the question of its accuracy). In PFA, 
one adds up "parallel plate" contributions to the force 
between two objects by including a force between each 
point on one surface and a corresponding point on the 
other surface, where corresponding points are connected 
by parallel "lines of interaction." In this geometry, we 
take the lines of interaction to lie parallel to the center- 
to-center displacement between two squares, as depicted 
by the inset in Fig. |3j but of course this choice is some- 
what arbitrary. (A similar choice was made by Ref. [26] 
to define the PFA force between two eccentric cylinders.) 
The PFA force between one pair of squares is then: 
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where D y / a 2 + (fr + s) 2 , H = 1 + s/h and A 
1 — s/\a\. The resulting net force is shown in Fig. |3j 
along with the contributions due to the isolated square- 
square and square-plate PFA forces (a separate line for 
the plate-plate contributions is not shown because this 
contribution is always very small). 

For comparison, Fig. [3] also shows the exact total force 
from Fig. [2j and it is clear that, while PFA captures the 
qualitative behavior of the oscillating force sign, in quan- 
titative terms it greatly overestimates the magnitude of 
the repulsive force. Of course, since it is an uncontrolled 
approximation in this regime there is no reason to expect 
quantitative accuracy, but the magnitude of the error il- 
lustrates how different the true Casimir force is from this 
simple estimate. The PFA estimate for the square-plate 
force, however, does help us to understand one feature of 
the exact result. If there were no plates, only squares, 
then the force would be zero by symmetry exactly at 
a = 0, and indeed the exact result including the plates 
has zero force at a « 0; clearly, the contribution to the 
force from the plates is negligible for a « 0, and this 
is echoed by the PFA P sq _ p i force. Also, using a PFA 
approximation, one can attempt to estimate the order 
of magnitude of the force contribution from the ends of 
the bracket, which was neglected in the exact calcula- 
tion. This contribution to the total force must decrease 
as ~ 1/L for a fixed a, and is estimated to be less than 
1% of the peak repulsive force for L > 60 A. 

Because the basic explanation for the sign changes in 
the force for this structure is fundamentally geometri- 
cal, we expect that the qualitative behavior will be ro- 
bust in the face of imperfect metals, surface roughness, 
and similar deviations from the ideal model here. The 



main challenge for an experimental realization (for ex- 
ample, to obtain a mechanical oscillator around the equi- 
librium point) would appear to be maintaining a close 
parallel separation of the brackets (although it may help 
that in at least one direction this parallelism is a stable 
equilibrium). Furthermore, although in this paper we 
demonstrated one realization of a geometry-based repul- 
sive Casimir force, this opens the possibility that future 
work will reveal similar phenomena in many other ge- 
ometries. 
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